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We study the dynamics of non-interacting quantum particles with two bands in the presence 
of random scattering. The two bands are associated with a chiral symmetry. After breaking the 
latter by a potential, we still find that the quantum dynamics is controlled by a non-Abelian chiral 
symmetry. The possibility of spontaneous symmetry breaking is analyzed within a self-consistent 
approach, and the instability of a symmetric solution is discussed. 

I. INTRODUCTION 

Recent experiments on two-dimensional two-band systems, such as graphene [J Q or the surface of 
topological insulators 01 , have revealed that transport in these systems is very robust with a universal 
behavior that does not depend on the sample quality. In general, the mechanism for the appearance of 
universal and robust properties in many-body systems is often associated with underlying symmetries 
and the spontaneous breaking of these symmetries. For instance, the long-range and long-time behavior is 
usually controlled by symmetries but it is not much affected by the fine details of the model. Apart from 
critical points in systems with second-order phase transitions, dynamical properties such as diffusion may 
also be related to symmetry properties. An example is a system of non-interacting Dirac fermions in two 
dimensions at the Dirac node with a random gap, where we have identified a non-Abelian chiral symmetry 
This symmetry transformation is usually called "supersymmetric" because it connects bosons and 
fermions [6| . It turned out that spontaneous breaking of this symmetry generates a massless fermion mode, 
caused by pairing a boson and a fermion. This is analogous to the Cooper pair in superconductivity which 
consists of a pair of fermions rather than a boson-fermion pair. Another difference is that the massless 
fermion decays as a power law and does not provide a coherent state. The physical effect of the massless 
mode is diffusion of fermionic quasiparticles. On the other hand, the system is insulating if the symmetry 
is not spontaneously broken. This can happen for a sufficiently large average gap Q. Moreover, diffusion 
and the metal-insulator transition are described by a simple scaling law which deviate from the Drude 
behavior of conventional metals Q . All this indicates that the spontaneous breaking of the non-Abelian 
chiral symmetry is the origin of the metallic behavior and the metal-insulator transition due to random 
gap opening of Dirac fermions. 

The idea of this work is to extend the analysis of the Dirac fermions to the more general case of tight- 
binding models with two bands. In particular, we will identify the non-Abelian chiral symmetry for a 
four-body Hamiltonian which is constructed from one-body Hamiltonians with broken chiral symmetry. 

II. MODEL 

In a chiral-invariant system states in the upper band at energy E are directly connected with states 
in the lower band at energy —E by a linear transformation U . In other words, the chiral Hamiltonian 
satisfies the anti-commutator relation 

[C/,i7]+-0. (1) 

For an eigenstate of energy E with H'^ e — E'^ e this relation implies that '^-e = U'^/e', i-e., the 
states of the negative (positive) spectrum are created by the linear transformation U from the states of 
the positive (negative) spectrum. A consequence of the vanishing anti-commutator is the chiral symmetry 
relation 

gaC/^gaC/ ^ (2) 



with a continuous parameter a. 
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A. Special Hamiltonians 



A two-band Hamiltonian can be represented in terms of Pauli matrices aj (j = 0, 3) as 

H = hoCTo + hiai + h2(J2 + ^^3^3 , (3) 

where hj are matrices in real space, and the associated eigenfunctions of H are spinor states. A symmetric 
example is a tight-binding Hamiltonian on a bipartite lattice, such as the honeycomb lattice. It can also be 
written in the form of Eq. ([3|), where the Pauli matrices refer to the two sublattices. For nearest-neighbor 
hopping on the honeycomb lattice the Fourier components of the coefficients read 

3 3 
hi = —t cos(Qj ■ k), h2 — —t sin(aj • k) 

with the lattice vectors oi ~ (— ■\/3/2, 1/2), 02 — (0, —1) and 03 — (■\/3/2, 1/2). ^13 either plays the role 
of a staggered potential that breaks the sublattice symmetry, or it represents a next nearest-neighbor 
hopping term on the honeycomb lattice 0. For /13 = the Hamiltonian has a chiral symmetry with 
U = while for /13 ^ the chiral symmetry is broken. We will treat Hamiltonians with broken chiral 
symmetry in Sect. IHI Al 

Remark: A prominent example for Eq. ([3]) is the two-dimensional Dirac Hamiltonian, where hj = idj 
(j — 1,2) and /13 = m is the mass term. It should be noticed that the Dirac Hamiltonian is not 
symmetric but chiral invariant for m = 0, and it satisfies Eq. ([1]) for [/ = 0-3 again, m ^ breaks the 
chiral symmetry. In that case the Dirac Hamiltonian satisfies 

aiH = -ffVi , (4) 

rather than the relation ([T]), with H* = , where with the matrix transposition Q. This relation is 
associated with particle-hole symmetry [loj . 



III. FEW-BODY APPROACH 



Now we consider briefly the dynamics of non-interacting particles in a random environment and derive 
a four-body Hamiltonian to describe transport behavior in the presence of random scattering. In order 
to study the motion of a quantum particle in space we calculate the transition probability of a particle 
going from site r at time t to site r' at time t' . The transition amplitude (or Green's function) for this 
process is 

Gr,t;r',t' = j e-^^(*-*')(F -E + ie)-l,dE , 
and the corresponding transition probability reads 

Pr,t;r',t' = \Gr,t;r' A'? = j g-'^^-^'^^*-*') {H - E + le);^, {H - E' ~ te);,\dEdE' . (5) 

Then the expansion of the wave function within the time interval [i, 0] can be written as 

, 2x ^ J2r'^^PrA;Q,Q 

A straightforward calculation for a translational invariant Hamiltonian H results in a quadratic increase 
of this expression with time t. This is ballistic transport, in contrast to diffusion, where it would increase 
linearly with time. Diffusive transport is possible when we consider a random Hamiltonian H and average 
the transition probability with respect to the distribution of H . 

The product of the two resolvents in Eq. ([5]) can be interpreted as a scattering process of two inde- 
pendent particles that move in opposite time directions because of ±ie (advanced and retarded Green's 
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functions). As a result, the dynamics is described by the scattering of two independent particles with iden- 
tical Hamiltonians. In principle, we could add more independent particles, since we always can project the 
resulting few-body Hamiltonian to the desired product of Green's functions. Using a functional-integral 
representation of non-interacting particles (cf. Appendix \^ , we can accommodate an arbitrary number 
of factors. It is important to notice that the average of the Green's function leads to an effective inter- 
action between the particles. We assume that the one-body Hamiltonians obey the relation ([IJ. Now we 
can take advantage of the relation ([T]) which allows us to switch between retarded and advanced Green's 
functions by applying the linear transformation U. This is discussed in the next Section. 



A. non-Abelian chiral symmetry of the four-body Hamiltonian 

We consider the extended Hamiltonian H + Hi, which is not chiral symmetric, and the transformation 
matrices Uj (j = 1, 2), for which we assume that they satisfy the (anti-)commutator relation 

[H,U,]+^[Hi,U,U^O (6) 

as a generalization of the relation ([!]). This includes the case Ui = U2 ^ U. Next we replicate the 
one-body Hamiltonian H + Hi to a two-body Hamiltonian H + Hi with 



H 



H 
H 



Hi 



Hi 

--Hi 



U 



Ui 
U2 



Together with Eq. ([6|) these expressions are connected by the relations 

[H,U]+ = [Hi,U]-=0- 
Now we introduce the four-body Hamiltonian 



H 



H + Hi 






H -Hi 



and the transformation matrix 



(7) 
(8) 

(9) 
(10) 

H . (11) 
Moreover, we notice that det{H + Hi-\-ie) = det{H — Hi+ie). Then we obtain for the graded determinant 




U = 



and obtain the anti-commutator relation 

This implies the non-Abelian chiral symmetry 

e^He^ 



ipiUi 






















. 





detg{H -I- ie) = detg(e^) = 1 . 
Eventually, there is a relation between advanced and retarded Green's function: 



U{H + iey^U-'^ 



i-H 



(H 



(12) 



Eqs. PU)) and (jlip are a non-Abelian generalization of Eqs. ([T]) and respectively. The symmetry 
transformation in Eq. ^ depends on the free parameters ipj , ip'^ . For the special case of a random 
Hamiltonian they are Grassmann variables (cf. Appendix El and Refs. [l,[l,[il|). This can be understood 
as if the upper block matrix (the lower block matrix) of the Hamiltonian ([5]) acts on bosons (fermions). 
Consequently, the off-diagonal blocks in the transformation matrix U transform bosons into fermions and 
vice versa. This requires that the parameters are Grassmann variables. 

Remark: The construction of the Hamiltonian H would work for H = 0. This implies that we could 
obtain the non-Abelian chiral symmetry even for a one-body Hamiltonian with only one band. The 
reason is that _ff is a two-body Hamiltonian with chiral symmetry. 
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IV. SPONTANEOUS SYMMETRY BREAKING 

The main question is whether or not the non-Abehan chiral symmetry in Eq. (jlip is spontaneously 
broken for the average Green's function {{H+ie}^^) after sending e — > 0. This can be studied by employing 
the self-consistent Born approximation (SCBA), which is equivalent to the saddle-point approximation of 
the functional integral in (|A3|) . For this purpose we introduce the (complex) self-energy 77 to approximate 
the average Green's function by 

{{H + z6)-i) « {{H) + ^e + tij)-^ = Go . (13) 
T] satisfies the self-consistent (or saddle-point) equation for disorder strength g 



where Xj are the Fourier components of hj of the Hamiltonian. The integration is performed with respect 
to the Brillouin zone, /i is a uniform chemical potential which replaces the Hi term in the Hamiltonian: 
Hi = /ido- This equation is equivalent to 

«Ai = 2 1 -1 + 5 / TTT TT— TTTT? ) ' (15) 



fe A? + A2 + A2 + z2 

which allows us to determine the two solutions z{±fi) for a given /x. The solutions of this equation can 
be distinguished according to their transformation behavior under the sign change fi — > — /i: There is a 
symmetric solution with 

z{~^J') = (16) 

which does not break the non-Abelian chiral symmetry and an asymmetric solution that obeys 

'^-'^ + ^ i Af + Al + k + zi-,r ) - {-'^' L A? + Ai 'xl + zi,Y ) ' 

where z{—iJLf' ^ The symmetric solution in implies z(0) = 0. Moreover, the effective 

symmetry breaking field is 

- ^ z{fi) + z{-^) ^ ry(/i) + ?7(-^) 

Whether the symmetric or the asymmetric solution is relevant depends on the stability of the solution. 
Inspection of the fluctuations around the symmetric solution, according to the expansion in Appendix |b1 
provides us a stable symmetric solution if the kernel of the quadratic form in Eq. (jBl|) is non-negative 
and unstable symmetric solution if it has negative eigenvalues. Thus spontaneous symmetry breaking 
is indicated by an instability of a symmetric solution. This will be discussed in the next Section for a 
specific choice of the Hamiltonian H. 

A. Example: linear spectrum 

In the case of a linear spectrum with s nodes, where we have Xj = kj {j = 1,2) with cut-off A (i.e. 
< fc^ -f fc| < A^) and A3 = 0, the integral in Eq. can be performed to yield 

r, = z^log(l + AV^') . (19) 

For /i = this equation has the symmetric solution z = ?; = and the asymmetric solution 

A 

Z ^ Tl = 

Ve2/'^9 - 1 ■ 
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For z 7^ we can rewrite (|19p as 

A2 = {^^hs^^'^i^ha-- _ 1^ . (20) 

Thus, the imaginary part of the right-hand side must vanish and its real part must be A^. There is 
no closed solution of these equations but for /i ~ we can expand the exponential function for the 
asymmetric solution and obtain a quadratic expression for the right-hand side 



Thus, the symmetry breaking field of ([T 




increases with /i. For small g we have /3/a ^ 1 such that 

(21) 

In physical terms, f] is the effective scattering rate in the Green's function (|13p . Our result means 
that the scattering rate increases as we go away from the nodes at /i = into the bands. This effect 
is not surprising, since the density of states also increases as we go away from the nodes. For this 
special case the eigenvalues of the stability matrix on large scales have been evaluated in Ref. Q as 
(l/g, 1/5 + 2, l/g — 2, 1/5). Thus, the symmetric solution is unstable for g > 1/2. 




V. SUMMARY AND CONCLUSIONS 



The starting point of our analysis is a one-body Hamiltonian H with chiral symmetry ([2|). After 
breaking this symmetry by a potential term, we have identified the non-Abelian chiral symmetry pil) 
of the corresponding four-body Hamiltonian ([5]). This four-body Hamiltonian can be used to describe 
random scattering processes, where the transformation from retarded to advanced Green's is performed 
by a chiral transformation. 

The non-Abelian chiral symmetry can be spontaneously broken after averaging the Green's function 
with respect to random fluctuations of the Hamiltonian. We have studied the general case in a self- 
consistent approximation and have identified symmetric and asymmetric solutions. For a special Hamil- 
tonian with linear spectrum around the nodes we have found that the symmetric solution becomes 
unstable for sufficiently strong random fluctuations. 

These results indicate that the symmetric two-band Hamiltonian has a similar non-Abelian structure 
as the 2D Dirac fermions with particle-hole symmetry. Like for the latter case, we anticipate that the 
spontaneously broken non-Abelian chiral symmetry leads to a massless fermion mode which describes the 
physics on large scales. 



Appendix A: Functional integral: distribution of the Green's function 



Now we combine the matrix structure of H in Eq. with the averaging procedure for a random 
Hamiltonian iJ. This will justify that the parameters in the transformation matrix U are Grass- 

mann variables. For this purpose we employ a functional-integral approach with bosons and fermions. 
The latter has been discussed extensively in the literature @ and specifically for the problem of two-band 
models in Refs. 0,i,[ni- Therefore, only a brief summary that is relevant for the symmetry discussion, 
is given here whereas details can be found in the mentioned literature. 
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We consider the Green's functions or products of Green's functions, averaged over a Gaussian distri- 
bution of some physical quantity inside the Hamihonian H. For simphcity we focus here on a random 
chemical potential and study {{H + je)~^), where we use the notation for the average 

{{H + le)-)) = 1: f{H + ^e)-l J] g-^'^/f d^, (Al) 

r 

with the normalization Af. Now we can apply a transformation from the random variable /i^ to the 8x8 
random matrix Qr 

with Hermitean 4x4 matrices Qr and Pr and 4x4 supermatrices 8, 8. Qr is related to the diagonal 
elements of the Green's function {H + ie)~l due to 

(iH + ie)-l) = / QrJlle-^'^a^^Q'^/^ViiQr] . (A3) 

Products of Green's functions can also be represented by this functional integral J the Jacobian 
of the transformation /Zr — > Qr and reads as a graded determinant 

J = detg((7?) +iefo + Q) . (A4) 

It is important to notice that J is invariant under the global symmetry transformation (|9]) because 

detg(e'^) = e^^^ = 1 with the graded trace Trg. Then we can apply the transformation Qr — e^QrC^ , 
which leaves the integral invariant in the limit e — >■ 0. By sending e to zero we are able to study 
spontaneous symmetry breaking within a saddle-point approximation of the functional integral (|A3p . 



Appendix B: Instability of the symmetric saddle point 

Fluctuations around the symmetric saddle point Q — can become unstable when there is an asymmet- 
ric saddle point. Gaussian fluctuations around the saddle-point solution indicate clearly the instability of 
the symmetric solution Q . The stability analysis can be applied to the integral in Eq. (|A3p by expanding 
the logarithm of the Jacobian in Eq. (jA4p around Q — Q: 



-Trg^iQl) 
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logdetg((i7) + ieio 



Q) = -Trg4{Ql) - -Trg^ Go-r.r'Qr'Ga-r' ,rQr 
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+ o(Q3). (B1) 



The symmetric solution is stable if the kernel of the quadratic form in Q is non-negative (i.e. it has no 
negative eigenvalues). On the other hand, it is unstable if one or more eigenvalues become negative. 
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